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Two-particle momentum correlations have been applied extensively to the determination of 
particle source distributions and flow effects in nuclear collisions [||]. HBT interferometry employs 
variations of the detected pair intensity in a two-particle momentum space to estimate the con- 
jugate space-time-momentum autocorrelation density of the emitting source. These established 
techniques depend directly on quantum-mechanical amplitude interference. We propose here an 
alternative source of momentum correlation induced by a combination of radial or Hubble flow and 
autocorrelation differences between particle pair types which should provide additional information 
on the two-particle space-time structure of the emitter. This new method is not primarily an in- 
terference effect. In the simplest interpretation it has a cosmological analog in the use of red shifts 
and a universal Hubble- flow hypothesis to map relative galactic positions. 

In this Letter we extend the usual treatment of two-particle momentum correlations to in- 
clude the possibility of non-chaotic or correlated particle emission from the hadronic freeze-out 
surface. We adopt a modified two-particle emission density which permits description of nontriv- 
ial differences between pair autocorrelations for different pair types. Pair position autocorrelation 
determines how nuclear flow is sampled by pair partners. If different pair types sample flow in sys- 
tematically different ways we predict a unique signature in two-particle momentum distributions. 

A general two-particle momentum distribution assuming independent pair emission can be 
written as P2(Pi,P2) = f d'^xid'^X2 S2{xi,X2,pi,P2) \^2{xi,X2,pi,P2)\^ , where |$2p may include a 
distorted- wave treatment to describe Coulomb and final-state interactions (FSI). In what follows 
we use X = {x\ + X2)/2 and y = xi — X2- If we assume a factorizable Wigner density (chaotic or 
independent particle emission), smoothness P, 0] and plane- wave propagation without final-state 
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interactions to a momentum detector we obtain 

P2(Pi,P2) ~ Piipi)Piip2)+ Jd*Xid^X2S,ix^,k)S^ix2,k)cos{q-y) (1) 

which is the conventional elementary expression for the sibling-pair (partners from same event) dis- 
tribution used in HBT analysis. The second term contains the Fourier transform of the configuration- 
space autocorrelation density. The condition / (fpi/Ei d^p2/E2 P2{pi, P2) = N{N — 1) can be used 
to normalize this distribution. 

We now seek a more general form of the two-particle Wigner density in which the momentum 
density is factorized as before but the emission density is not factorized as usual. We assume instead 
that pair emission varies locally with partner separation distance in a way which may depend on 
some property of the pair, such as total isospin for example. One can think of this as the effect of 
an initial-state interaction rather than a final-state interaction. Guided by FSI treatments it is 
then natural to factor the two-particle emission density g{xi,X2,pi,P2) into dependence on mean 
pair position x and on partner separation y, retaining the possibility of parametric dependence on 
momenta. We then obtain 

S2{xi,X2,pi,P2) ~ g+ix,k)g^{y,q) fi{xi,pi) fi{x2,P2)- (2) 

The argument of fi{xi,pi) /i(x2,P2) = Ei E2 exp ^_ P'^M^i)+P2-u{x2) ^ g rearranged as pi ■ 

u{xi) + p2 ■ u{x2) = k ■ + -^q ■ Au, where we adopt the definitions k = {pi + P2)/2 , q = Pi — P2 
and define Sn = u{xi) + u{x2), Au = u{xi) — u{x2). Assuming a longitudinally-invariant radial 
flow velocity field given near z = by u(x) = jix) (1, /3_|_(x)) Q and particle transverse momentum 
near Y = given hy p = {mt,Pt) we have p ■ u{x) = 7(x) (mt — P±{x) ■ pt). Further assuming 
that 7(2;) ~ 1 over most of the volume we can approximate the sum and difference terms as 
k-T,u = niti + 10112 - k • {I3^{xi) + /3_l(x2)) and g ■ = -q • (/3_l(xi) - (3^{x2)). 

For discussion purposes we represent P2(pi) P2) as a sum of two terms ^(k, q) -|- -B(k, q), where 
the second (interference) term contains the usual HBT Fourier transform (= in the case of mixed 
pairs). The first term, which we now examine in detail, is ordinarily written as the product of two 
single-particle distributions. We now approximate radial flow observed in HI collisions Q by analogy 
with cosmological Hubble flow and assume that V/3 = H X from which follows A I3{x) ~ H Ax, 
that is, relative transverse velocity is proportional to transverse separation for any two points in 
the emitting system. This approximation may be unrealistic at the scale of a nuclear diameter 
but quite reasonable at the hadron scale of primary interest here. Combining results we have for 
^(k, q) the product of three factors 

^(k,q) = mtimt2e t ■ d xg+{x, k) exp i—k ■ xj ■ d y g^{y,q) exp i— y j{3) 

The first two factors are equivalent to the single-particle radial-flow case for the particle pair as a 
unit. The third factor contains potentially new information about the pair autocorrelation density 

9-{y,q)- 
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The two-particle correlator ||5|, |^ can be expressed as C2(k, q) = a(k, q) -|- 6(k, q), where for 
instance a(k, q) = ylo;,j(k, q)/Are/(k, q) and 6(k, q) contains the HBT Fourier transform from 
which the large-scale source geometry is inferred. In the ratio a(k, q) the first two factors of Eq. 
@ cancel, leaving the ratio of two integrals on particle separation. The standard assumption of 
chaotic particle emission implies that the autocorrelation densities for object and reference systems 
are equal {gobj{y,Q) ~ 9refiy,Q)) and therefore a(k, q) « 1. We argue that this assumption is 
arbitrary and that the consequences of gobj 7^ Qref should be tested experimentally. For the ratio 
we then obtain (we have changed notation: gobj) 

ref 

a(k,q) ^ JdS{9obj{y,q) - gref{y,q)}cosh(^^q-yj (4) 

where V is the volume inferred from the autocorrelation distribution, and we argue that the sink 
term in the integrand is suppressed by symmetry, leaving cosh as the dominant term. 

Put simply, just as radial or Hubble flow causes broadening (blue shift) of the single-particle pt 
distribution and two-particle k distribution to higher momenta it also distorts the distribution on 
pair momentum difference q. However, in the latter case the degree of distortion depends on the 
autocorrelation density on partner separation as opposed to mean position on the emitter. If this 
autocorrelation is different for object and reference distributions (Hubble flow is sampled differently 
in the two cases) the difference is reflected to lowest order in a large-scale quadratic departure from 
unity in the ratio of the two distributions on q. Again, this is not primarily an interference effect. 

We can model the difference gobjiy, q)—gref{y, q) by assuming that event-mixed (reference) pairs 
follow an autocorrelation density on the transverse emitting space described by a gaussian with rms 
transverse radius R reflecting the size of the emitting system. We further assume that for sibling 
pairs the autocorrelation density may be modified by a second smaller concentric gaussian with 
radius r < R, added or subtracted (correlation or anticorrelation) depending on pair properties, 
and with amplitude e <C 1 reflecting the strength of this pair (anti) correlation. Normalization of 
densities is maintained. Source radii inferred from HBT analysis, sometimes referred to as 'lengths 
of homogeneity' , more generally correspond to the second moments (covariance matrix) of the 
autocorrelation density which may have an arbitrary shape near the origin, thus precluding use of 
the curvature matrix |^] to represent source sizes. 

The leading term after expanding the cosh in Eq. (^) contains the factor < (q-y)^ >Ag- We can 
express the dot product in terms of longitudinal and transverse components for y and q referred 
to mean vectors x and k respectively. The cross term linear in qi yi is zero by symmetry. This 
leaves the two terms < (qtyt)'^ > and < {qiyif' >■ We assume that the autocorrelation difference 
/S.g{y, q) localizes yt and yi to a common rms value y ^ r, and that /^g{y, q) may contain a strong 
correlation between (j)y and 4>q, the opening angles between position and momentum vector pairs 
respectively, in which case we obtain < (q • y)^ >Ag~ r^{4A;^ < sin? {(p/2) > W < cos^{(p/2) >}, 
where k and q are now the scalar mean and difference oi pti and pt2, and (p is the azimuthal opening 
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angle between momentum vectors. Gathering results we obtain 



a{ptuPt2) ~ l^e^{r/Rf{l-{r/Rf}{2k'{l-<cos{^)>^g) 

+ (1+ < cos(0) >Ac,)/2} < cos{(P) >Age [0, 1] (5) 



where P± is the rms flow speed for the gref{y,Q) distribution, and < cos{4>) >Ag represents the 
autocorrelation-averaged pair momentum opening angle. The negative resp. positive sign corre- 
sponds to correlation resp. anticorrelation of sibling pairs on separation distance relative to a mixed 
reference. Thus, any difference in the autocorrelation density of emitted particles for sibling pairs 
relative to a mixed reference (breakdown of chaotic-source assumption) as measured by a difference 
in rms partner separation between sibling and mixed pairs combined with local Hubble flow should 
result in a contribution to the correlator which is quadratic in pt mean and/or difference. Eq. (|5|) 
should model any emission scenario from uncorrelated emission throughout a transparent volume 
to normal particle emission from an opaque cylinder. 

The value of <cos{4>) >Ag could be established from experiment by comparing the amplitudes 
of observed quadratic dependences on k and q. Due to the finite mean free path of pions in 



nuclear matter and possible opacity effects |1C] the mean value of the cosine may be significantly 
different from zero, and possibly close to unity. For example, the relationship between partner 
separation and opening angle for normal emission from an opaque cylinder of radius R is cos{<j)){y) = 
cos{2cos~^ y^l — {y/2R)'-^}. This expression remains close to unity (> 0.9) even for substantial 
separation {y/R < 0.5). The implications of opacity for source space-time structure is a very active 
area of research. Information on opacity has been derived by comparing Rout and Rgide obtained 
from HBT analysis. However, such comparisons are ambiguous due to the well-known cross talk 
between source duration and thickness in the 'out' direction. Determination of <cos{(f)) >Ag by this 
new method could provide important independent information on the issue of opacity. Currently 
there is unresolved conflict between models describing a 'sudden transition to full transparency' 
1^] (< cos{(/)) >Ag~ 0) and models describing sustained emission from a thin surface layer on an 



opaque source |10|] {<cos{4>) >Ag~ !)• Alternatively, the value of <cos{(p) >Ag may be influenced 
by resonance decay kinematics. 

In summary, the large-f; region of two-particle momentum space may carry significant dynamical 
information, may in fact be oversubscribed by different dynamical phenomena at some level of 
sensitivity. Analysis of the two-particle momentum correlator in this region could provide detailed 
information on the emitting surface beyond the large-scale features of the space-time-momentum 
density revealed by established HBT interference techniques. We here propose a new strategy 
combining nuclear Hubble flow with a generalized emitter autocorrelation density to probe details 
of the hadronic freeze-out surface at small length scales, especially to explore the dependence 
of the emission process on particle-pair characteristics such as total isospin. We find that any 
departure from the chaotic emission assumption used in standard HBT analyses combined with 
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local Hubble flow should result in a quadratic deviation from unity of the two-particle correlator at 
large transverse momentum sum and/or difference. By observing such correlations one may learn 
about two-particle spatial and angular correlations at small length scale on the emitter which could 
indicate whether hadrons are generated from a bulk prehadronic state (QGP) or decouple with 
correlations typical of a hadronic cascade. 
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